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Abstract. A special Infeld-HuU factorization is given for the Askey- Wilson 
second order g-difference operator. It is then shown how to deduce a gener- 
alization of the corresponding Askey- Wilson polynomials. 
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1. Introduction. 

A given second order g-difFerence operator 

H{z; n) = U2{qz; ra)E^ + ui{qz; n)Eq + uo{qz; n) (1) 

where 

^\{'^{z;n)) = '^{zq';n), zeR,i,neZ, (2) 

is said to be factorizable according to Infeld-HuU (this is a g-version) if and 
only if the following product can be performed [7], [15]: 

H{z;n) - fi{n) = {ri{z;n)Eq + g{z;n)) {C{z;n)Eq + f{z;n)) 
H{z;n+l)- ii{n) = (C(^;n)E, + /(^;n)) (77(^;n)E,+5(^;n)) (3) 

The operators ri{z] n)'Fjq+g{z] n) and C,{z] n)'Fjq+f{z] n) are called " lowering" 
and "raising" respectively. It is thus a particular case of the well known 
Darboux transformation (where n is no more a variable but all simply an 

"Permanent address : Universite du Burundi, Faculte des sciences, Departement de 
Mathematique, B. P. 2700 Bujumbura, Burundi, East Central Africa. 



1 



index) [5], [14], [18], [21]. 

Supposing that U2 doesn't depends on n, let us reduce (3) in the following: 

H{z;n) - fi{n) = {u2{z)Eg + g{z;n)) {u2{z)Eg + f{z;n)) 
H{z;n+1) - ll{n) = {u2{z)Eg + f{z;n)) {u2{z)Eg + g{z;n)), (4) 

where H{z;n) = U2{z)H {z; n) . The corresponding system then reads: 

f{qz;n+l)+g{z;n+l) = f{z-n)+g{qz-n) 

f{z-n+l)g{z-n+l) = f{z;n)g{z;n) + fi{n)-il{n+l) (5) 

As a general property, the system (5) (more precisely, a system equivalent 
to (5)) is known to be a discrete version of the Toda lattice (see for exam- 
ple [23]). It is easy to see that the intertwining relations in (4) (isospectral 
deformations) will then constitute a discrete version of the corresponding 
Lax-pair. The system (5) can also be considered as a discrete delay-Painleve 
I [6]. It is morever known to admit symmetries [20]. We are however in this 
work concerned in applications of the IHF rather than in study of its intrin- 
sic structure. 

An alternative to the IHF is the Inui factorization [8] where the "up- 
stair" (raising) and "down-stair" (lowering) operators are not obtained from 
a direct factorization of the studied second order difference operator itself 
but from some difference or contiguous relations satisfied by its eigenfunc- 
tions given in a general a priori known form. An Inui type factorization for 
the Askey- Wilson operator [1] (see [10] for the canonical form) can be found 
in [16], [9] and [4] where interesting applications were also obtained. In [22] 
and [23] IHF techniques were applied to the three-term recurrence relations 
for orthogonal polynomials. 

In [19] and [3] were given special IHF of the following difference operator 

CT(s)AV + f(s;/)A- A(n;/) (6) 

solutions of which are the /-th difference of the discrete hypergeometric poly- 
nomials on a linear lattice Pi'\s) := A'(P„(s)) [17]. Here A = T - 1, 
V = 1 — (T'p(s) := p{s + i)) , (T and f being polynomials (in s) of 
degree < 2 and 1 respectively, A a constant (in s). In [19] the role of the 
"variable of factorization" (as n in (3)) is played by the order / of the cor- 
responding difference (ra remaining fixed) while in [3] that role is played by 
the degree of the corresponding polynomials (/ remaining fixed (to zero, for 
simplicity)). Evidently, the IHF from [19] can be extended to the case of 
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non-linear lattices. In that case, one intersects (formally) the evoked above 
Inui type factorizations from [16], [9], [4]. 

The main objective of the present work is to give an extension of [3] to the 
case of non-linear lattices. But here, to be concrete, we will work with the 
Askey- Wilson second order g-difference operator in its canonical form [10], 
instead of the Nikiforov-Suslov-Uvarov hypergeometric difference operator 
[17]. A simple and nonsignificant reduction of the lattice and an elementary 
change of variables lead indeed at the latter to the former. The IHF that 
we are going to carry out for the Askey- Wilson operator will furnish various 
formulas such as strict (i.e. without any perturbation of parameters) dif- 
ference relations, recurrence relations and Rodrigues type formulas for the 
general sequences of eigenfunctions of the Askey- Wilson operator. Next, a 
somewhat converse reasoning will suffice to attain the most ambitious (see 
the conclusion) aim of this work: to implement the IHF system (5) as an 
independent source of the Askey- Wilson Polynomials ^. All those discus- 
sions will be held in the following section. In the last one, we will extract a 
Rodrigues type formula allowing to produce (starting at the difference hy- 
pergeometric functions obtained in [2]) a sequence of functions generalizing 
the Askey- Wilson polynomials. 



2. The factorization method for the Askey- Wilson second or- 
der g-difference operator. 



z-\-z 



The Askey- Wilson polynomials (AWP) [1], Vn{x{^))} x{^) — 2 
fined by: 



are de- 



„ . . ^^ {ab,ac,ad;q)n , ( q '"^ahcdq" ^,az,az ^ 

Vnixiz)) - h 



q;q 



(7) 



where the basic hypergeometric (or g-hypergeometric) series rft^s read: 



ai, . . . , Or 



bi,...,b,s 



q-z 



00 / 



k=0 



(5i, . ..,bs;q)k 



^^!iik(_i)(i+«-oy 1+^-0(2) 



{q;q)k 



^The original idea of generating discrete orthogonal polynomials from the IHF system 
(5) comes from [15] but the approach adopted there (as in [7]) already doesn't allow to 
handle totally the Hahn case. 
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with 

{ai, . . .,a.r;q)k ■= {ai; q)k ■ ■ ■ {a-r; q)k, 

while 



k-i 

{a; q)o := 1, {a; q)k := H " ^^'')' ^ = 1, 2, 3, . . . (8) 

8 = 

They satisfy the following second order g-difference equation [10]: 

£Vn{x{z)) = Kn)Vn{x{z)) (9) 

where 

£{z) = v{z)E, - {v{z) + v{z-')) + v{z-')E-' (10) 

with 

(1 - a^)(l - bz){l - cz){l - dz) 
""^^^ ~ {I - z^){l - qz^) 

E;(^n(x(^))) = VMq'z)),ie Z, (11) 

is the Askey- Wilson second order g-difference operator 
and 

A(n) = -(1 - g-")(l - ahcdq^-^). 

In anticipation on the coming g-integration and basing ourself on the ex- 
perience from [3], we are led to reformulate, before the factorization, the 
equation (9) as follows: 

{A{qz)El - [A{qz) + B{qz) + }C{qz)\{n)]E, + B{qz)}Vn{x{z)) = (12) 

where 

A(^) — A-2Z-^+A-iz-^+Ao+Aiz+A2Z^ 
^ ' qz—z~^ 

N _ A2Z-^+Aiz-^+Ao+A-iz+A-2Z^ 

jC{z) = z- z-^ 

A-2 = 1; A-i = — {a -\- b -\- c -\- d); Aq = ab -\- ac -\- ad -\- be -\- bd -\- cd 
Ai = — {abc + abd + bed + acd) ; A2 = abed 

Consider next the following g-difference operator 

n{z; n) = A{z) [A{qz)El - [A{qz) + B{qz) + IC{qz)X{n)]Eq + B{qz)] . (13) 
The central result of this work is the following 
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Proposition 1 The operator 7i in (13) factorizes into 

%{z; n) - ii{n) = (yl(^)E, + g{z; n))(yl(^)E, + T{z; n)) 
n{z; n + 1) - i^{n) = {A{z)Eg + T{z; n)){A{z)Eg + g{z; n)) (14) 

whith T and Q such that the class of the Askey-Wilson polynomials becomes 
invariant under the action of the "raising" and "lowering" operators. 

Proof. The operatorial relations (14) are equivalent to the system: 

T{qz;n) + g{z;n) = -{A{qz) + B{qz) + IC{qz)X{n)) 
T {z; n)Q [z; n) = A{z)B{qz) — ii{n) 
^q{T{z-n)-g{z-n)) = (A(n + 1) - A(n))/C(g^) (15) 

where Aqh{z) = h{qz) — h{z). 

Using the g-integration, we first transform the system (15) in : 



:F{z; n) + :F{qz; n) = -{A{qz) + B{qz)) - (/3_i - AM)^-i _ f]^ 

-{Pi + X{n)q)z 

T {z; n)Q [z; n) = A{z)B{qz) — ii{n) 

Q{z-n)-T{z-n)= Y:^\^rz' (16) 

where /3_i = MliiLLizMli) ^ = gP-i', remaining arbitrary for the moment. 
Observing the first equation in (16) and then using the last one, it becomes 
sensible to search J-{z; n) and then Q{z] n) under the forms: 

jrt^.^y^ F-2z-''+F_iz-^+Fo+Fiz+F2z' 



g(^^.^y^ (F_2-P-i)z-^ + (F-i-Po)z-^+Fo + (Fi+poq)z+(F2+piq)z^ l^-^rj^ 

Taking ^u, /3o, -f-2j -f^-ij -foj -f^ij -^2 as unknowns, the system (15) will then 
be transformed in an algebraic system of 16 equations for 7 unknowns. To 
solve it (by hand), one first determines those unknowns from 7 equations 
and then very delicately ensures himself that they satisfy all the remaining 
9 equations. 
The result is: 

F-2{n): Mn)-jMn+i) _ ^ (^g) 



F2(n) : —r-i 



A(»)g-A(»+1) 2 _ q^+qA2 
y g + 1 
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F-i{n) 

Fi{n) 
Fo{n) 
li{n) 



r/o A(n)g-A(n+l) 2 , A(n+1)-A(n) 2 

(A(n)-A(n+l))g3tl^ l_g2 t/ i_g t/ 



2 



Ai+qA_l 



9 + 9" 



g/?o("-) 

2 2 

?3 - Ao(g + g2) + A2(g - 1) + g2(A(n) + A(n + 1))} 



Aq + + AoA^q-^ + Fo(n)/3_i(n) + F_i(n)/3o(n) 

-2F_2(n)Fo(n)-F2^(n). 



Thus, (17) is determined and a solution of the IHF problem (14) is found. 
It remains to prove the invariance of the class of the AWP under the action 
of the "raising" and "lowering" operators thus obtained. From (14), the 
intertwining relations follow: 



'H{z-n + l){A{z)'E.q + T{z-n)) 
'H{z-n){A{z)'Eq + g{z-n)) 



{A{z)'E.q + T{z-n))'H{z-n) 
(yl(^)E, + g{z- n))n{z; n+1) (19) 



Let us note that those relations, as well as the IHF (14) are in fact valid 
Vra G Z. So, if for a given number no £ Z , a given function {^) is an 
eigenvector of the operator 7i{z;no), corresponding to a given eigenvalue v 
(take z/ = 0, in our situation), then the sequence of functions 



k-i 



no +k ( 



n {A{z)Eg + T{z- no + i))*„o {z),k= 1, 2, 3, 



(20) 



8 = 



are respectively eigenvectors of 7i{z; no + k) corresponding to the same eigen- 
value z/ = 0. The following formulas will be useful: 



T{z; n) - Q{z; n - 1) = 

[2A2g-lg" - 2q-^]x{z) - Mn)+Mn-1) ^ (y^ ^ 



{A{z)E, + :F{z; 0))(1) = [2A2q-' - 2]x(^ 



/3o(0) 



+ 



-qA. 



2q 



(21) 

(22) 



Considering (20) and (14), we obtain the following difference relations (we 
let tHr^J^) ■.= ^no{z)y. 

^*n+l(^) 

-//(n)^*„(^) 



{A{z)E, + T{z;n))^^n{z) 
{A{z)Eg + g{z;n))^^n+i{z) 
ra = no, no + 1, no + 2, . . . 



(23) 
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Now, using (23) and (21), we obtain the following three-term recurrence 
relations: 



Considering the particular case of (20) when {^) is the AWP of order 
no, Vnoixi^)) 1 the invariance of the class of the AWP under the "raising" 
operation means that ^no+k{z) will be a non vanishing constant multiple of 
Vno+k{x{z)), the AWP of order no + k, k>l. 

It follows from (22) and (24) (with no = 0) that there exists a sequence of 
polynomial (in xi^)) eigenfunctions (of degrees 0,1,2,3,...) of the Askey- 
Wilson operator £, invariant under the "raising" operation. On the other 
side, according to (23), the same sequence is invariant , up to a multipli- 
cation by a constant, under the "lowering" operation (as long as ra > 0). 
But as the constant multiple of AWP are the unique non trivial polynomial 
eigenfunctions of the Askey- Wilson operator, this suffices to show the cited 
invariance and the proposition is completely proved. 

Thus, if in (20) we consider the case of AWP and set no = 0, we obtain 
the Rodrigues type formula: 



c{n)Vn{x{z)) = l[[A{z)Eg+T{z;i)]{l),n=l,2,3,... (25) 



c{n) being some no null constant (in z). It is not difficult to see that (25) is 
equivalent to 




(24) 



n = no + l,no + 2, . . . 



n-1 



1 



n-1 



d{n)Vn{x{z)) 




+ .F(^;^)](/>(^)),n= 1,2,3,... 



(26) 



where 



A{z) 



a formula similar to that from [3]. 
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Now, remembering that the relations (19) as well as the factorization 
(14) are valid Vra G Z ^ let us give the following notes for completness: 
Remark first that as in (20), the functions 
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^n,-x-h[z) := n('^(^)Eg + a(^; no - l-i))*„o(^),A; = 0,1,2,... (27) 



8 = 



are respectively eigenfunctions of the operators 7i{z;no — l — k). They satisfy 
the difference relations ( letting ^"^noi^) •= ^^noi^))'- 

-fi{n)^^n+i{z) = {A{z)Eq+T{z;n))^^n{z) 

^^n{z) = {A{z)Eg + g{z;n))^^n+i{z) (28) 
n = no — 1, no — 2, . . . 

and the recurrence relations 

fi{n)^'<i'n+i{z) +^ *„-l(^) = 

([■2,-" - 2A,_rVM^) + -" )'>!•..(.) (29) 

n = no — 1, no — 2, . . . 

For rao = and ipo{z) = 1 (polynomial case), (27) gives an "extension" of 
(25): 



^^.(n+i){z) := n('^(^)Eg + 0{z; -(^ + 1)))(1), n = 0, 1, 2, . . . (30) 

8 = 

But a direct verification of the relation 



{A{z)E, + g{z;-l)){l) = (31) 

shows, without surprise, that the functions produced in (30), "extending" 
the AWP relatively to n (as an index, from to Z~), are all null (~''^'_(„_|_i) (z) = 
0, ra = 0,l,2,...). 

The possibility of unifying the "raising" and "lowering" operations in (25) 
and (30) in types (23) or (28) pair of difference relations (and so in types 
(24) or (29) recurrence relations) with n = . . . — 2, — 1, 0, 1, 2, . . ., is subor- 
dinated to the constraint /«( — 1) = and this is actually the case. 
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We are now led to make a somewhat converse reasoning: 
We remark first tliat 



T{z;-l) 

Q{z;-i) 







B{qz) 

A{z) 



(32) 



Tliis leads together with the first equation in (15) to a new form of the 
equation (12): 



{g{qz- -1)E2 - [T{qz- n) + g{z- n)]E, + T{z- -l)}Vn{x{z)) = 0. (33) 



We recall that this equation is verified by T and Q given by (18) and (17) 
and the AWP given by (25) or (26) for ra > 1 while Vq = 1 (the unit is a 
solution of (33) for ra = thanks to the first equation in the system (5)). 
Let us return now the situation taking the IHF system (5) as a starting 
point. Our point resides in the following 

Proposition 2 Consider the IHF system (5) with the unique initial condi- 
tion /u( — 1) = and let T and Q he any of its solutions, then the equation 



{g{qz; -1)E2 - [T{qz; ra) + g{z; ra)]E, + T{z; -l)}Vn{z) = (34) 



admits a sequence of solutions given by a "Rodrigues type formula" (similar 
to (20)) applied to Vq = 1 and satisfying srict difference relations and three- 
term recurrence relations. Morever, if we add the condition 



for some constants (in z) cq and ci, the obtained sequence is of polynomial 
(in x{z)) type. 

Proof. To prove this propostion, one needs essentially to remark that if T 
and Q are solutions of the system (5) with initial condition /*( — 1) = 0, then 
the operator 



n{z; ra) = g{z; -1) [g{qz; -l)El - [:F{qz; ra) + g{z; n)]E, + :F{z; -1)] (36) 



admits a type (14) IHF. The remaining part of the reasoning is similar to 
that of the preceding proposition, reason for which we omit it and we assume 
that the proposition is proved. 



J^{z; ra) - g{z; ra - 1) = co(ra)x(2:) + ci(ra) 



(35) 



9 



It is clear that if we add the analogous of the two first initial conditions 
in (32), we obtain the AWP. In other words, we remark in our satisfaction 
that the IHF system (5) can be considered as an independent source of the 
AWP. 

3. The generalization of the Askey- Wilson polynomials. 

We envisage now to show how to produce a sequence of functions ^n{z,t) 
generalizing the AWP in that sense that they are eigenfunctions of the 
Askey- Wilson operator £ (see (10)) and ^n{z,t) — > T^nixi^)) while t — > 1. 
We showed above that if the eigenvalue of the Askey- Wilson operator is 
A(ra) = — (1 — g~")(l — a5c(ig"~^), then an IHF as in (14) is realizable. Now, 
using (15), one finds that conversely, the possibility of a type (14) IHF im- 
plies necessarily that the eigenvalue in action is *A(ra) = — (1 — ig~")(l — 
abcdt~^ q"^~^) . On the other side, it is clear that the presence of the parame- 
ter t in *A(ra) does not in no way annoy the IHF (14) (thanks to the specific 
dependence of *A(ra) at t). This means that if in the expressions (18), we 
substitute A(ra) by *A(ra), the resulting expressions for T and Q, say and 
solve a type (14) IHF problem for an operator, say ^7i{z;n), obtained 
from 7i{z;n) replacing A(ra) by *A(ra). 
Let us consider next the equation: 

£y{z,t)=' X{0)y{z,t), (37) 

or equivalently 

'n{z-0)y{z,t) = 0. (38) 

This equation is not new: A little wider equation (where the Askey- Wilson 
operator is replaced by the corresponding Nikiforov-Suslov-Uvarov operator 
[17]) has been explicitely solved in [2]. In our situation, according to [2], one 
first solves the non-homogeneous equations: 

\{0)^ u{z,t,a) = G{z,t,a) (39) 

where 

{at,abcdat-^q-^;q)c^ _i 

G{z,t,a) = —— —{az,az ;g)oo (40) 

a(aba, aca, ada, aq; q)oo 

and ((t; g)oo := limj_!.oo (f; g)i, for the following values of a: 

q q q 

a := 1 a := — a := — a := — -. 

ab ac ad 
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The corresponding solutions are the functions [2]: 



(az,az ,q)oo y-^ (at, abcdat q ,aaz,aaz ,q)i ; 

u{z,t,a] = - — - > — ; q. (41) 

[aaz,aaz ,q)^~^ [aba, aca, ada, aq; q)i 

Let us remark that the AWP appear here as the functions u{z,q~"',l), 
ra = 0,l,2,... 

The solutions of the homogeneous equation (37) or (38) are then obtained 
by opering adequate (considering (40)) linear combinations of any two of 
the functions (41). 

Among the solutions of the homogenous equation (37) or (38), we are inter- 
ested in those having a constant limit when t converges to the unit. Let us 
take the "adequate" linear combination of u{z,t, 1) and u{z,t,a), a / 1: 

{t, abcdt~^q~^ , aba, aca, ada, aq; q)oo 
y{z, t) := u{z, t, 1) - a abcdat-^q-^- q) ^' " ^^^^ 

It is easily seen that y{z, 1) = 1. However, if y{z,t) is another "adequate" 
linear combination of u{z,t,ai) and u{z,t,a2) where ai and a2 differ both 
from the unit and each other, then y{z, 1) / constant. In other words y{z, 1) 
is the nonconstant solution of the Askey- Wilson equation for ra = 0. 
Our point here resides in the following 

Proposition 3 If y{z,t) is the solution in (42), of (37) or (38), then the 
functions 

<^o{z,t) := y{z,t), 

n-l 

^n{z,t) := Y{[A{z)^g+' T{z;i)]y{z,t), (43) 

8 = 

n= 1,2,3,..., 

generalize the AWP in that sense that they are eigenfunctions of the Askey- 
Wilson operator £ (corresponding to *A(ra) ) and they converge to them when 
t converges to the unit. 

Proof. As already noted, the operator admits a type (14) IHF but 
now with T and Q replaced by and ^Q. Consequentely, the fact that 
the functions in (43) are eigenfunctions of the Askey- Wilson operator £ 
corresponding to *A(ra), is a consequence (thanks to type (19) intertwining 
relations) of that ^o{z;t) is an eigenvector of £ corresponding to *A(0). To 



11 



be assured that the functions in (43) converge to the Askey- Wilson polyno- 
mials (of course, up to a multiplication by a constant), we need essentially 
to remember that *A(ra) converges to A(ra) when t converges to the unit and 
then compare the r.h.s. of (43) and (25), which proves the propositon. 

For completness. It's perhaps worth noting that an iterative applica- 
tion of the "lowering" operator A{z)'Eq -|-* Q{z; 1)), i = 0, 1, 2, ... to 
y{z,t) leads to an evident "extension" relatively to n (as an index, from 
to Z~) of the functions <i>„(2:,i) in (43), generalizing thus the already 
evoked identically null functions (30). Morever if ^^'^(2:) (or ^^n{z)) are the 
eigenfunctions obtained by applying ra-times the "raising" (or "lowering") 
operator to y{z, 1), so their generalizations are obtained by performing sim- 
ilar operations starting now at y{z,t). Let us note finally that the present 
generalizations are also difference hypergeometric functions (in the sense of 
[2]) reason for which they are closely related to the AWP. 

Conclusion. We have (mainly) shown explicitely in this work and in [3] 
that all the classical (up to the AWP) orthogonal polynomials can be recov- 
ered solving the IHF system (5) (ra being presumed to be the degree of the 
searched polynomials). We have serious indications allowing to believe that 
there (i.e. a little modification of the IHF system (5)) is exactely the central 
"front door" when attempting to handle concretely non-classical classes of 
orthogonal polynomials such as the discrete semi-classical orthogonal poly- 
nomials or even the discrete Laguerre-Hahn orthogonal polynomials (the 
largest nowadays known class of orthogonal polynomials) [11], [12], [13]. We 
reserve a concrete discussion in this matter for a later publication. 
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